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ABSTRACT 


An  infrared  technique  has  been  used  to  non  i  tor  the  the  ml  decomposition  of 
carbonyl  fluoride  (GJF2)  in  argon  and  nitrogen  diluent  behind  incident  and  n- 
flected  shock  waves.  Data  was  taken  in  the  temperature  range  2400°  to  36000  * 
at  total  pressures  between  0.2  and  26  atn.  Direct  sanpling  of  reflected  shock 
nixtures  with  a  ti*e-of- flight  mass  spectrometer  provided  knowledge  of  the  de¬ 
composition  products.  The  pressure  and  tenperature  dependence  of  the  CDF2 
dissociation  rate  constants  are  discussed  in  terns  of  the  Rice-Raasperger- 
SLassel  (Iff)  uninolecular  theory.  Rate  constants  for  the  fluorine  extraction  reac¬ 
tions  and  CJF  disproportionation  reaction  were  obtained  by  cu’rve-fitting  the 
complete  COF2  kinetic  histories  with  computed  profiles.  A  chemical  nonequi librium 
stream-tube  program  which  includes  wall  boundary  layer  effects  was  used  for  data 
analysis. 
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I. 


INTRODUCTION 


At  present  the  chemical  rates  for  fluorocarbon  decomposition  and  oxidation 
kinetics  are  needed  for  input  to  a  number  of  aerospace  engineering  problems. 

Hass  spectrometer  study  on  the  pyrolysis  of  Teflon  plastic  (polytetrafluoroethv- 
lene)  in  atmospheres  of  air  and  oxygen*  shows  that  carbonyl  fluoride  (COF2)  is  an 
important  low- temperature  (800°  K)  oxidation  intermediate  and  that  at  higher 
temperatures  carbon  dioxide  (CO2)  and  perfluorome thane  (CF4)  dominate  the  reaction 
products.  Measurements  to  identify  the  infrared  and  visible  spectrum  of  an  cir- 
Teflon  boundary  layer  in  a  hot  subsonic  arc  jet  have  shown  the  presence  of  the 
CN-violet  and  red  bands  and  structure  in  the  2  to  6p  region  characteristic  of  the 
COF2,  CO  and  CO2  molecular  bands. 2*3  Single-pulse  shock  tube  experiments  of  the 
tetraf luoroethylene  oxidation  in  the  range  of  1200°  to  2000°  K  show  product  distri¬ 
butions  of  C0F2»  CO,  CF4,  C2F5,  and  small  amounts  of  C3F6  and  C02.4  In  a  shock- 
tube  study  of  the  difluoromethylene-oxygen  reaction,  infrared  measurements  have 
directly  identified  COF2  as  the  initial  oxidation  product  of  these  reactants. 
Preliminary  results  show  that  at  temperatures  <2200°  K,  the  COF2  builds  up  to  a 
1  steady  concentration  behind  the  shock  front  and  tnat  above  2700°  K,  a  nonequili- 

V  brium  "overshoot"  in  COF2  takes  place  in  the  reaction  mixture. 5 


From  a  measurement  of  the  thermal  decomposition  rates  of  molecules,  the  rate  con¬ 
stants  for  their  formation  can  be  obtained  by  the  principle  of  detailed  balancing, 
i.e.,  taking  the  ratio  of  the  dissociation  rate  to  the  equilibrium  constant  of  the 
reaction.  This  approach  has  previously  been  used,  for  example,  to  determine 
fluorine-fluorocarbon  dissociation  and  recombination  rate  constants  of  CF4  and 
CF3.^  In  the  piesent  report  another  part  of  fluorocarbon  chemical  kinetics  has 
been  investigated  and  deals  with  the  thermal  decomposition  of  carbonyl  fluoride. 

A  shock  tube  is  employed  to  dissociate  C0F2  in  dilute  mixtures  of  argon  or  nitro¬ 
gen.  The  rate  of  COF2  decomposition  is  measured  by  monitoring  the  COF2  fundamental 
band  at  5.25 ft  with  a  filtered  InSb  infrared  detector.  Rate  constants  for  COF2 
initial  dissociation  are  obtained  over  a  pressure  range  from  about  0.2  to  26  atm 
between  2400°  and  3600°  K.  The  pressure  and  temperature  dependences  of  these 
rate  constants  are  expressed  in  terms  of  the  RRK  theory  for  unimolecular  reactions. 
In  addition,  the  complete  COF2  decomposition  records  are  analyzed  with  a  chemical, 
nonequilibrium  shock-tube  program  to  obtain  rate  constants  for  the  fluorine  ex¬ 
traction  and  recombination  reaction  of  the  C0F2  decomposition  mechanism. 
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II.  EXPERIMENTAL  PROCEDURE 


A.  SHOCK  TUBE 

The  COF2  decomposition  experiments  were  conducted  behind  incident  and  reflected 
shock  waves  in  an  emission  spectroscopy  shock  tube  and  behind  reflected  shock 
waves  in  a  shock  tube  coupled  to  a  time-of-flight  mass  spectrometer. 

The  optical  shock  tube  was  constructed  of  1.5-inch  l.D.  stainless  steel  tubing 
and  was  designed  with  a  3. 5- foot  driver  section  and  a  10-foot  driven  section. 

The  shock  tube  and  associated  gas  handling  system  were  pumped  down  to  less  than 
3  x  10”  3  torr  with  a  water-cooled  oil  diffusion  pump;  they  had  a  leak  rate  of 
about  10”  3  torr  per  minute.  Cold  driving  with  hydrogen  of  helium  against  Mylar 
diaphragms  of  2-mil  to  5-mil  thickness  was  used  to  generate  shock  waves  into  the 
test  chamber.  The  shock  velocity  was  measured  with  a  series  of  platinum  heat- 
transfer  gages  situated  along  the  driven  section.  The  output  of  these  gages  were 
displayed  on  a  time-mark  folded  oscilloscope  sweep,  operated  by  a  Radionics  (Model 
TWN-2A)  triangular  wave  and  marker  timing  generator.  Transit  times  between  stations 
were  measured  to  within  +1  psec.  The  maximum  error  in  the  shock  velocity  was  about 
+0.3  percent  and  introduced  an  uncertainty  of  less  than  +15°  K  in  the  calculated 
incident  shock  temperature  or  +40°  K  for  the  reflected  shock.  Infrared  emission 
from  the  shocked  gas  was  monitored  with  a  filtered  InSb  detector  (Type  ISC-301) 
and  was  viewed  through  calcium  fluoride  windows  perpendicular  to  the  shock-tube 
axis. 

The  mass  spectrometer  shock-tube  apparatus  was  constructed  from  1-inch  l.D.  Pyrex 
glass  pipe  and  has  a  5-foot  driver  section  and  a  15-foot  driven  section.  A  Bendix 
(Model  14-206)  time-of-flight  mass  spectrometer  is  attached  to  the  shock  tube 
through  the  "fast  reaction  chamber."  The  reflected  shock-heated  gas  was  sampled 
at  the  shock-tube  endwall,  which  is  a  flat  stainless  steel  plate  having  in  its 
center  a  hyperbolic  nozzle  with  a  2-mil  inlet  and  10-mil  exit  (Englehard  Indus¬ 
tries,  Englehard,  New  Jersey).  Mass  spectral  recordings  were  taken  at  25-//sec 
intervals  just  before  and  after  shock  arrival.  The  time-resolved  ion  peaks  were 
displayed  on  a  Tektronix  535A  oscilloscope  in  combination  with  a  type  CA-pre- 
amplifier  and  multiscan  generator. 7  For  initial  shock-tube  pressures  of  3  to  8 
torr,  the  background  pressure  in  the  mass  spectrometer  before  shock  arrival  ranged 
between  10”6  and  10“5  torr.  A  complete  description  of  the  mass  spectrometer  shock- 
tube  apparatus  may  be  found  in  an  earlier  report.® 

B.  REAGENTS 

Infrared  measurements  of  COF2  decomposition  were  made  behind  incident  argon  and 
nitrogen  shocks  containing  COF2  in  a  mole  ratio  of  2:100  with  the  inert  gas.  The 
same  measurements  were  made  behind  reflected  shocks  using  0.25:100  and  0.5:100 
COF2  argon  mixtures.  Mass  spectral  data  on  the  COF2  decomposition  reaction  were 
taken  with  a  5:100  C0F2-argon  mixture.  The  COF2  gas  used  in  this  study  was  ob¬ 
tained  from  Peninsular  Chemresearch  (Gainesville,  Florida)  and  had  a  stated  purity 
of  95  percent.  Infrared  analysis  of  COF2  completely  decomposed  in  excess  argon 
showed  that  the  commercial  gas  contained  about  4.4  +2.4  percent  CO2.  The  buffer 
gases  were  of  research  grade  (argon,  99.999  percent  purity  and  N2,  99.997  percent 
purity) . 
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III. 


INFRARED  MEASUREMENTS 


The  Infrared  absorption  spectra  of  the  COF2  molecule  is  found  to  consist  of  a 
number  of  vibrational-rotational  bands  between  2m  and  38m. (9)  Recently  in  this 
laboratory,  spectral  band  intensities  in  the  2p  to  6fi  region  of  CUf^  (Figure  1) 
have  been  determined  by  employing  a  "fast  scanning"  monochromator  to  record  the 
infrared  spectra  of  reflected  shock  heated  samples  of  COF2  in  excess  argon. (10) 

For  the  purpose  of  studying  the  thermal  decomposition  kinetics  of  COF2,  the  funda¬ 
mental  vibrational  band  at  5.25^  was  monitored  with  a  filtered  InSb  detector  to 
observe  the  time  rate  of  change  of  the  C0F2  concentration  behind  the  shock  wave. 
Examples  of  Infrared  oscillogram  records  depicting  C0F2  decomposition  at  several 
shock  temperature  conditions  are  shown  in  Figure  2.  The  COF2  concentrations  of 
the  study  were  typically  less  than  2  x  10“?  mole/cc  for  which  the  infrared  emission 
at  5.25 m  was  below  8  percent  blackbody  intensity  over  the  temperature  range  of  the 
experiments.  Under  these  conditions,  the  C0F2  gas  was  considered  optically  thin 
and  the  measured  infrared  signal  was  taken  to  be  linearly  dependent  on  concen¬ 
tration.  From  the  infrared  emission  peak-height  at  the  shock  front  and  the  initial 
concentration  of  COF2  behind  the  shock  wave,  each  oscillogram  record  was  self¬ 
calibrating.  It  was  found  in  the  experiments  that  the  CO  decomposition  product 
and  the  trace  amount  of  CO2  impurity  in  the  COF2  reaction  mixture  radiated  in  the 
bandpass  of  the  5.25m  filter.  Corrections  to  the  COF2  infrared  records  were  made 
by  subtracting  out  the  CO  (less  than  10  percent)  and  C02  (less  than  5  percent) 
spectra.  Necessary  correction  factors  were  obtained  from  measurements  of  the  in¬ 
frared  signal  generated  by  known  quantities  of  CO  and  CO2  in  argon  shocks. 
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Figure  1  AB0LUTE  INTENSITY  OF  THE  C0F2  EMISSION  SPECTRA 
IN  THE  2fi  TO  REGION 


EMISSION 


TEMP  2I85°K  PRESSURE  1.07  atm 
50  fi  sec  BETWEEN  VERTICALS 


TEMP  2460°K  PRESSURE  1.23  atm 
50  fi  sec 


TEMP  3040  °K  PRESSURE  l.56atm 
20  /x  sec 


TEMP  3250 °K  PRESSURE  1.60  atm 
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88-7011 

Figure  2  EXAMPLES  OF  COF?  INFRARED  EMISSION  AT  5.25„  BEHIND 
INCIDENT  SHOCK  WAVES  INTO  2:100  COF2-ARGON  MIXTURES 
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IV. 


HASS  SPECTROMETER  MEASUREMENTS 


Reaction  Mixtures  of  reflected  shock-heated  COF2  in  excess  argon  were  sampled 
directly  through  the  shock-tube  endvall  and  analyzed  in  a  time-of-f light  mass 
spectrometer.  The  ness  spectral  data  were  used  to  study  the  change  In  C0F2  mass 
peak  with  tine  (25-psec  intervals  after  shock  compression)  and  to  determine  the 
nature  of  the  decomposition  products  to  gain  insight  about  the  C0F2  decomposition 
mechanism.  A  typical  mass  spectra  of  a  COF2  decomposition  mixture  is  seen  in 
Figure  3.  Ion  peaks  for  COF2C66),  COF(47),  C02(44),  Ar(40),  Ar(36),  and  CO(28) 
are  clearly  identified.  Interestingly,  mass  peaks  for  CF3169),  CF2(50),  and 
CF(31)  are  absent  in  the  spectra.  The  histories  of  the  COF2  ion  peak  were  used 
directly  to  obtain  the  initial  dissociation  rate  constants  for  the  COF2  molecule. 
The  mass  spectrometer  records  show  that  CD  is  the  major  molecular  end-product  of 
COF2  decomposition.  This  observation  is  in  agreement  with  CD  infrared  analysis 
at  4. Cp.  It  is  found  that  after  complete  decomposition  of  COF2,  the  CD  concen¬ 
tration  in  the  mixtures  equals  the  stoichiometric  CO  content  of  the  COF 2  molecules 
present  at  the  beginning  of  the  reaction. 


T 
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V.  THERMAL  DISSOCIATION  RATE  CONSTANT  OF  CARBONYL  FLUORIDE 


The  kinetics  of  COF?  thermal  dissociation  were  studied  over  a  pressure  range  from 
about  0.2  to  26  ata  in  argon  shocks,  and  froa  about  0.4  to  1.5  ata  in  shock-heated 
nitrogen.  The  teaperature  interval  of  the  data  extended  froa  2400P  to  3600°  K. 

The  rate  constant  for  thermal  dissociation  of  C0F2 


Ar,N2 

R1  C0F2  — ►  COF  +  F 

was  determined  froa  the  initial  slope  of  the  infrared  eaission  records,  corres¬ 
ponding  to  the  first  10-percent  decay  in  the  COF2  concentration.  Sone  dissoci¬ 
ation  rate  constants  were  also  oeasured  off  Bass  spectrometer  oscillogram  records 
for  50  percent  COF2  dissociation.  The  "apparent”  dissociation  rate  constant  de¬ 
fined  by 


dlcOF,! 

U'OF,ldt  k°b' 


(1) 


where  (  COF25 is  the  COF2  concentration  behind  the  shock  wave,  is  plotted  against 
reciprocal  teaperature  in  Figures  4  and  5.  The  average  total  concentration  associ¬ 
ated  with  each  set  of  dissociation  rate  constants  is  denoted  by  a  normalized  den¬ 
sity  equal  to  the  density  of  the  shocked  gas  divided  by  the  density  of  the  gas  at 
1  ata  pressure  and  273°  K.  Examples  of  the  shock  tube  experiments  performed  in 
this  study  are  given  in  Table  I. 


It  is  seen  froa  the  data  of  Figures  4  and  5,  that  the  apparent  first-order  rate 
constants  are  neither  independent  nor  linearly  dependent  on  total  concentration, 
but  rather  shew  a  systematic  fall-off  with  concentration  at  a  given  teaperature. 
This  type  of  behavior  for  thernal  dissociation  reactions  of  polvatomic  molecules 
has  been  explained  by  a  number  of  theoretical  formulations.  (11-13)  The  COF2 
dissociation  rate  constants  of  this  work  are  fitted  by  the  general  rate  constant 
equation  derived  froa  the  Rice-Ramsperger-Kassel  (RRK)  theory  on  unimolecular  re¬ 
actions.  The  general  rate  constant  formula  of  the  RRK  theory  is  given  by 


(2) 


RT 

where  *  Ac  is  the  Arrhenius  high-pressure  rate  constant.  Ec  is  the  minimum 
energy  for  dissociation.  A  is  the  frequency  factor.  The  quantities  m,  c,  and  Z, 
are  the  number  of  effective  oscillators  participating  in  the  dissociation  process, 
the  total  concentration,  and  the  effective  collision  frequency,  x  is  the  dimension¬ 
less  variable,  (E-E0)'RT.  The  RRK  rate  constant  for  C0F2  decomposition  com¬ 
pared  with  the  general  rate  constant  from  Slater's  theory. 
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TABLE  I 


r 


cof2  dissociation  in  argon  shocks 


Experiment 

No. 

Shock 

Temperature 

<°K) 

Total 

Concentration 

(oole/cc) 

Apparent 

Rate  Constant 
(sec*1) 

Total 

Pressure 

(ata) 

2:100  COF2  “  Argon 

1* 

2405 

6.245  x  lO-6 

1.11  x  103 

1.23 

2 

2545 

6.270  X  10~6 

2.42  x  103 

1.31 

3 

2880 

6.345  x  20~6 

9.45  x  103 

1.50 

4 

\ 

3195 

6.395  x  10~6 

2.42  x  104 

1.68 

\ 

i 

i 

0 

.5:100  C0F,  -  Argon 

i  5** 

2678 

1.586  x  10~5 

4.94  x  103 

3.48 

6 

2880 

2.231  x  10^ 

1.29  x  104 

5.27 

i 

;  7 

j 

3304 

1.686  x  10~5 

6.79  x  104 

4.57 

1 

;  8 

3680 

1.300  x  10~5 

1.24  x  105 

3.93 

9 

i 

j  2475 

3.840  x  10~5 

2.23  x  10? 

7.80 

| 

-5 

3 

1 

10 

2655 

3.930  x  10  ^ 

6.70  x  10 

8.60  ; 

-5 

u 

11 

3010 

4.090  x  10 

3.62  x  10 

10.10  \ 

12 

3270 

4.180  x  10~5 

1.08  x  105 

11.2C  j 

0.25:100  C0F2  -  Argon 


!  13 

2343 

1.19  x  10"4 

1.62  x  103 

- — i 

22.9  ! 

-4 

3 

14 

2410 

1.20  x  10 

1.93  x  10 

23.8  j 

\ 

15 

2795 

1.18  x  10~4 

2.71  x  104 

27.2  j 

-4 

4 

i 

16 

2905 

1.12  r.  10 

4.91  x  10 

26.7 

17 

3100 

1.10  x  10~4 

9.88  x  104 

28.0 

♦Incident  shock  data 
♦♦Reflected  shock  data 


TABLE  I  (Cood'd) 


COF2  DISSOCIATION  IN  NITROGEN  SHOCKS 
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where  #  is  related  to  the  mean  frequency  and  asplitude  factors  of  vibrational 
nodes  leading  to  reaction.  The  paraneters  of  Equation  (2)  giving  a  best  fit  to 
the  observed  COF2  rate  constant  were  evaluated  as  follows.  The  icv  *s  were  deter- 
cined  by  plotting  the  experinen tally  observed  rate  constants  at  a  given  temperature 
against  pressure  (Figure  6)  and  to  a  first  approxiaat ion  least  squares  fitting  the 

p 

data  with  the  phenoiaenological  equation,  kot>s  - - ,  based  on  the  Lindeaann 

p  •  t 

aechaniss  for  a  uniaolecular  reaction* The  k^  fs  were  next  used  to  deteroine 
A  by  a  least  squares  fit  of  the  Arrhenius  rate  constant,  taking  EQ  =  91.912  kcal/ 
nole  (see  Figure  7).  The  RRK  rate  integral  Equation  (2)  was  solved  on  an  IBM 
360  computer  for  A  «  6.14  x  1011  sec"1  and  Ec  *  91.912  kcal/aole  where  ■  and  Z 
were  varied  systematically  until  a  satisfactory  fit  was  achieved  with  one  set  of 
the  experinental  data  (Figure  8).  The  solid  curves  in  Figures  4  and  5  are  the 
theoretical  RRK  rate  constants  and  show  that  with  the  kinetic  paraneters 


A  Ec  Z 

Gas  sec"1  kcal  nole  a  cc /carle  sec 

Argon  6.14 (1011)  91.912  6  7.5  1011 

X2  6.14 (1011)  91.912  6  1.0  1011 

an  accurate  correlation  (within  40  percent)  of  both  the  argon  and  nitrogen  data 
is  obtained  over  the  concentration  range  of  the  study.  However,  it  is  pointed 
out  here  that  these  values  of  m  and  Z  are  not  at  all  unique  in  fitting  the  data. 
Reasonable  curve  fits  within  the  scatter  of  the  experimental  points  arc  also  ob¬ 
tained  for  m  -  4  and  Z  «  5.0  x  1011  cc/oole  sec.  Calculations  with  the  RRK  inte¬ 
gral  form  -6  equally  can  show  agreement  with  the  experimental  data. 


•In  order  to  be  consistent  with  current  JAXAF  thennochemical  data,  the  first  and  second  bond  dissociation  energies  of 
COFs  were  based  oc  the  enthalpy  change  for  the  reaction  CO2  =  CO  -  2F,  -VlT  =  160.'  kcal  mole  and  were  taken  to  be 

in  the  same  ratio  as  the  first  and  second  bond  dissociation  energies  oi  CF^.  Analogous  to  P(F  —  CF^l  D(F  —  (  F,  >  - 
I.'  *,  the  values  oi  l'(F  -  <  OF!  and  [>fF  —  (  O'  ate  91 -rl 2  and  6S.1'  kcal  mole,  respectively. 
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PRESSURE  (atm) 

M-9929 


Vertical  bar  indicates  experimental  scatter  of  data.  Horizontal  bar  shows  pressure 
interval  of  data  point.  Number  of  experiments  used  to  evaluate  data  point  near  a 
given  pressure  appears  beside  vertical  bar.  Solid  tines  are  nonlinear  least 
squares  fit  of  data  with  k^  =  k^  p/p  +  r. 


Figure  6  PRESSURE  DEPENDENCE  OF  COF2  DISSOCIATION  RATE  CONSTANT 
WITH  TEMPERATURES 
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•tc 


Vertical  bar  shows  statistical  spread  of  k_  values. 


Figure  7  LEAST  SQUARES  FIT  OF  EXTRAPOLATED  LIMITING  HIGH-PRESSURE 
RATE  CONSTANT  WITH  ARRHENIUS  FORMULA 
FOR  E0  =  91912  CAL/MOLE 
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VI.  MECHANISM  AND  KINETICS  OF  CARBONYL  FLUORIDE  DECOMPOSITION 


Mass  spectrometer  and  infrared  evidence  show  that  under  the  experimental  con¬ 
ditions  of  this  study  COF2  would  mainly  decompose  to  CO  and  fluorine  atoms.  A 
set  of  reactions  which  would  lead  to  these  results  are  the  thermal  dissociation 
reactions  of  COF2  and  COF,  the  fluorine  extraction  reactions,  and  the  COF  dis¬ 
proportionation  reaction  (Table  II).  The  rate  constant  of  R1  was  obtained,  as 
described  in  the  preceding  section,  from  the  initial  decay  in  the  infrared  records 
of  COF2  behind  the  shock  front.  The  thermal  dissociation  rate  constant  of  the 
COF  radical  (R2)  is  based  on  the  difference  in  the  first  and  second  bond  dissoci¬ 
ation  energies,  D(F-COF)  =  91912  cal/mole  and  D(F-CO)  *  68570  cal/mole  of  the  COF2 
molecule.  The  rate  constants  for  the  fluorine  extraction  reactions  and  COF  dis¬ 
proportionation  reaction  were  derived  from  an  Arrhenius  expression  by  incorporat¬ 
ing  the  energy  of  reaction  at  0°  K  in  the  exponential  term  and  curve  fitting 
measured  COF2  profiles  over  the  complete  chemical  relaxation  period  to  establish 
the  magnitude  of  the  pre-exponential  factor.  Calculated  COF2  profiles  are  ob¬ 
tained  from  a  nonequilibrium  shock-tube  computer  program  which  solves  simultaneously 
the  Rankine-Hugoniot  differential  equations  for  a  normal  shock  wave  and  the  appro¬ 
priate  chemical  kinetic  equations.  The  present  shock-tube  computer  program  also 
includes  the  effect  of  wall  boundary  layer  on  the  temperature  and  density  behind 
the  shock  front  (Appendix  a).  Species  profiles  typical  of  a  chemically  relaxing 
shock  wave  into  a  2:100  C0F2-Ar  mixture  are  displayed  in  Figure  9.  Also  depicted 
in  the  figure  are  the  chemical  rates  of  the  decomposition  reactions.  It  is  seen 
that  for  these  particular  experimental  conditions,  the  COF2  thermal  dissociation 
rate  (Rl)  dominates  the  first  20  percent  of  COF2  decomposition  and  then  later  R5' 
becomes  the  most  important  path  for  COF2  removal.  The  reasonableness  with  which 
the  chemistry  and  rate  constants  in  Table  II  can  predict  the  experimental  COF2 
decomposition  profiles  for  a  number  of  experimental  shock  conditions  is  demon¬ 
strated  by  the  results  in  Figure  10.  The  solid  curves  in  the  figure  are  computed 
with  the  non-equilibrium  shock-tube  program  including  wall  boundary  layer  effects. 
The  broken  curve  is  calculated  without  the  boundary  layer  corrections  to  temper¬ 
ature  and  density  behind  the  shock  front.  For  the  experimental  conditions  chosen, 
the  boundary  layer  effect  is  negligible  near  the  shock  front,  but  does  have  a 
noticeable  influence  on  the  analysis  of  the  data  taken  at  long  times  behind  the 
shock  wave.  In  Figure  11,  various  chemical  models  for  COF2  decomposition  are  com¬ 
pared  to  assess  their  compatibility  with  the  experimental  data.  According  to  the 
sensitivity  of  the  rate  data  and  the  results  shown  in  Figure  11,  the  rate  con¬ 
stants  of  Table  II  ar^  considered  to  be  the  most  probable  values  within  a  spread 
of  a  factor  of  ten. 


SPECIES  CONCENTRATION,  moteyfcc 


LABORATORY  TIME,  microseconds  LABORATORY  TIME,  microseconds 

88-7012 

#,  experimental  C0F2  infrared  measurements.  Solid  lines  were 
calculated  with  chemistry  and  rate  constants  of  Table  II. 


Figure  9  flONEQUILIBRIUM  STREAM-TUBE  SPECIES  PROFILES  BEHIND 
INCIDENT  SHOCKWAVE  IN  A  2:100  COF2-ARGON  MIXTURE 
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♦  INITIAL  TEMP  -  END  TEMP  (®K) 
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VII. 


DISCUSSION 


Apparent  rate  constants  for  COF2  dissociation  between  2400°  to  3600°  K.  exhibited 
a  narked  fall-off  as  the  pressure  of  inert  gas  was  increased  over  a  100-fold 
range.  The  Arrhenius  pre-exponential  factor  A  obtained  by  extrapolating  the  cx- 
perinental  data  to  the  high-pressure  uninolecular  limit  appears  to  be  in  agree¬ 
ment  with  values  (10*2  to  10**  *ec~A)  usually  considered  "normal**  for  uninolecular 
reactions.  The  nunber  of  "effective"  oscillators  used  in  the  MtlC  rate  integral 
to  fit  the  apparent  first-order  rate  constants  was  chosen  to  be  6,  which  is  the 
total  nunber  of  normal  nodes  in  CDF2  and  the  maximum  nunber  permitted  by  theory 
for  this  type  of  molecule.  Values  of  a  less  than  6  were  not  considered  because 
they  result  in  fits  whose  temperature  dependences  are  qualitatively  steeper  than 
the  trend  indicated  by  the  data.  Use  of  values  greater  than  6  would  have  been 
contradictory  to  theory.  Although  the  kinetic  parameters  expressing  the  experi¬ 
mental  dissociation  rate  constants  of  this  work  my  be  somewhat  vagarious,  non- 
theless,  it  will  be  interesting  to  use  then  as  exanples  to  treat  other  theoretical 
relations  which  might  want  to  be  compared  with  experiment. 


2  /8*kT\* 

For  hard  spheres,  the  classical  collision  frequency  Z  is  Oj?  - j  which  Is 

9.1  x  1012  cc/mole  sec  for  COF2 -argon  partners  and  1.1  x  10*3  for  COF2-N2  partners 
(the  collision  diameters  are  taken  to  be  aAr»  3.4A,  *  3.71A  and  «C0F2  *  4.6A, 

the  nean  diameter  of  CO7  and  CF*).^1-*^  comparison  of  the  hard-spheres  collision 
frequency  with  experiment  yields  steric  factors  of  0.008  and  0.009  for  COF2 
dissociation  by  argon  and  nitrogen.  The  steric  factor  generally  is  associated 
with  the  fraction  of  collisions  that  cause  dissociation  of  critically  energized 
nolecules  and  is  considered  to  be  a  measure  of  the  steric  hindrance  in  the  col¬ 
liding  pair. 

Theoretically  the  transition  pressure  for  a  uninolecular  reaction  is  approximately 
related  to  the  high-pressure  frequency  factor  A  by 

A[\E  (E0  *  AE)]“  RT  =  PZ  (4) 

where  at  moderate  pressures  AE  on  the  average  has  the  value  (m  +  1)  RT.^^  With 
E0  «  91.912  kcal/mole  and  the  experimental  values  of  m  -  6,  A  -  6.14  x  10AA  sec”A, 
the  transition  pressure  for  the  uninolecular  dissociation  of  COF2  in  argon  at 
3000°  K  is  calculated  to  be  15.4  atm,  which  appears  to  be  in  line  with  the  ob¬ 
served  pressure  dependence  (Figure  6). 

The  concept  that  in  a  unimolecular  reaction  the  activated  complex  is  in  quasi¬ 
equilibrium  with  the  normal  molecule  has  led  to  the  formulation  of  Eyring's  tran¬ 
sition  state  theory  and  a  rationale  for  the  Arrhenius  pre-exponential  factor. 

The  quantity  A  is  related 
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where  Q*  is  the  partition  of  the  activated  coapla  and  Q  is  the  parti  tier;  function 
of  the  nonzal  molecule.  Theoretically  to  evaluate  the  pre-exponential  factor  for 
the  OOF 2  unlmolecular  dissociation,  the  electronic  contribution  to  the  partition 
function  for  both  the  complex  and  normal  molecule  is  assumed  to  be  unity.  The 
transition-state  of  the  activated  COF2  molecule  is  next  thought  of  being  a  OCF 
fragment  bounded  by  a  weakly  attached  fluorine  atom.  It  is  also  assumed  that  the 
OOF  fragment  rotates  and  vibrates  independently  of  its  satellite  fluorine  atom. 

On  the  basis  of  the  model  presented,  the  rotational  partition  function  of  the 
activated  complex  was  calculated  with  the  moments  of  inertia  for  a  CDF  radical. 

The  ratio  of  the  rotational  partition  function  of  the  activated  couples  to  that 
of  the  normal  molecule  was  thus  calculated  to  be  0.348.  Moments  of  inertia  for 
COF2  nod  COF  were  obtained  from  JAXAF  the  rmoc  hem  leal  tables17.  In  its  activated 
state,  the  deformation  modes  of  the  COF 2  molecule  are  assumed  to  be  drastically 
weakened  in  which  case  the  vibration  partition  function  of  the  complex  is  largely 
determined  by  the  V},  V2,  V4  modes  of  the  COF 2  molecule,  that  is,  the  in-phase 
(C-F),  the  COO)  vibration  and  the  out-of-phase  (C-F)  vibration.  Accordingly, 
the  ratio  of  the  vibrational  partition  functions  was  calculated  to  be  0.316.  Sub¬ 
stitution  of  these  numbers  for  the  ratio  of  the  partition  functions  in  Equation  (>) 
gives  a  calculated  pre-exponential  factor  of  6.86  x  1011  sec”1  at  3000°  K.  Although 
the  result  is  in  remarkable  agreement  with  the  experimental  value,  the  coincidence 
is  believed  to  be  more  fortuitous  than  real,  since  the  nature  of  the  complex,  though 
reasonably  assumed,  must  be  regarded  as  only  a  guess. 

It  is  concluded  from  mass  spectra  obtained  in  the  present  study  that  the  fluoro¬ 
carbon  products  CF4,  CF3,  and  CF2*  which  are  produced  by  C2F4  oxidation,  do  not 
appear  to  form  by  COF 2  decomposition.  Infrared  measurements  of  4.3a  radiation 
also  did  not  show  kinetic  evidence  for  CO2  production,  nor  did  ultraviolet  measure¬ 
ments  at  2660  A  indicate  the  existence  of  CF2  species  in  the  COF 2  decomposition 
mixtures.  In  accordance  with  these  observations,  the  reactions 

COF,  +  COF,  — CO,  +  CF, 

2  2  2  4 

C0F2  +  COF  —  C02  +  CF3 

COF  -f  COF  — ►  C02  +  CF2 

cof2  +  m  — ►cf2  +  o  +  m 

would  not  appear  to  be  important  in  the  chemistry  for  C0F2  decomposition. 
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ATPEXDIX  A 


FOMXATIOS  OF  A  GENERALIZED  CHEMICAL  WKEQll  LIBRIUM  STREAM 
TOE  COMPUTER  PROGRAM  WITH  IttLL  ROUCO*  T  LATER  EFFECTS 


Recently.  Teare"  has  pointed  oat  that  in  shock-tube  chemical  kinetic  studies,  the 
accuracy  of  the  results  nay  be  significantly  impaired,  if  variations  in  the  flow 
*  vity  and  state  properties  due  to  wall  boundary  layer  growth  are  neglected  in 
->  ta  analysis.**  Generally  it  is  understood  that  the  presence  of  wall 
.  ry  layer  In  the  chock  tube  causes  the  gas  behind  the  shock  wave  to  behave 
■*  if  it  were  passing  through  a  divergent  channel  where  the  flow  velocity  de- 
i-r<- :  and  the  gas  density  and  thermodynamic  temperature  increases  with  increas- 

iTjf  distance.  An  analytical  basis  for  these  effects  is  provided  by  the  theoretical 
work  of  Mirels.  who  has  node  use  of  Duff's  concept  of  a  ""limiting  separation 
distance"  to  obtain  a  local  similarity  solution  for  the  boundary  layer  displace¬ 
ment  thickness.19  The  constraints  placed  on  the  problem  require  that  the  mass 
per  unit  tine  integrated  over  the  entire  shock  tube  area  be  conserved  and  that 
the  mass  if  gas  entering  the  shock  front  be  equal  to  the  mass  entering  the 
k'WJdary  layer  at  the  limiting  separation  distance. 

Fc:  the  purpose  of  analyzing  shock-tube  chemical  kinetic  data  taken  at  long  tines 
behind  the  shock  front,  a  cnnequi librium  stream-tube  computer  program  has  been 
developed  to  include  area  variation  due  to  wall  boundary  layer  growth.  The  quan¬ 
tities  describing  the  flow  in  wave-stationary  coordinates  are  shown  in  Figure  A-l 
which  is  similar  to  the  diagram  given  in  Reference  19.  The  "jump"  conditions 
across  the  shock  wave  are  obtained  from  the  Hugoniot  conservation  equations, 

nass:  p,L  t  = 


t  2 

nfcraentun:  Pj  -  Pj  L*s*  "  Peo  *  ?co  uro~ 


•>  2 
enthalpy:  hj  -  1  -  LV  =  heo  '  1  -  uco 


The  equations  governing  the  chemical  relaxation  zone  behind  the  shock  wave  are 
expressed  in  time  differential  form  by 


continuity: 


d(puA ) 

^T 


P'uc  a 

p\  Us  A0 


0 
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A  SUBSONIC  flow  cross-sectional  area, 
pre-exponential  factor 

c  SPED ES  CONCENTRATION 

E  activation  energy 

h  STATIC  ENTHALPY 
N  GRAM  MOLECULAR  WEIGHT 
«■  TEMPERATURE  EXPONENT 
P  STATIC  PRESSURE 
R  IDEAL-GAS  CONSTANT 
T  THERMODYNAMIC  TEMPERATURE 
t  TIME 

U  SUBSONIC  FLOW  VELOCITY 


1  CONDITION  AHEAD  OF  SMOCK  WAVE 
0  SHOCK  FRONT  CONDITION 
•  CONDITION  BEHIND  SHOCK  WAVE 
f  FORWARD  REACTION 
•>  PROPERTY  OF  »*TH  SPECIES.  ETC 
i  i’TH  REACTION 
t  LABORATORY  COORDINATE 
P  PARTICLE  COORDINATE 
r  REVERSE  REACTION 


U,  SHOCK  WAVE  VELOCITY 
«  DISTANCE  BEHIND  SHOCK  WAVE 


.f  LIMITING  SEPARATION  DISTANCE 
p  DENSITY 

»•  REACTANT  STOICHIOMETRIC-COEFFICIENT 

Figur.  A— 1  SHOCK  TUBE  FLOW  DIAGRAM  FOR  FiXED  SHOCK 
WAVE  COORDINATE  SYSTEM  INCLUDING  WALL 
BOUNDARY  LAYER  DEVELOPMENT 
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For  a  laadnar  boundary  layer  (Reynolds  maber  less  chan  3  x  10?)  the  area  varia- 
tlon  of  the  divergent  channel  fro*  Mirels*  results  is  given  by 

Afl  'A  =  1  -  (x'*f)°-5 
The  ideal-gas  equation. 

p  *  ('£ci/2>)RT- 


is  used  to  relate  the  thermodynamic  temperature  to  pressure  and  density.  The 
molar  enthalpy  of  each  species  is  expressed  in  terns  of  temperature  by  a  power 
series. 

hi  -  [LlT  .  ~  L2T?  *  j  t,T»  -  ^L4T<  -  L,T-‘  -  L6]i  . 

in  vhicl«  the  temperature  coefficients  are  determined  from  curve  fits  of  JAXAF 
thermochemical  data.  The  variation  in  species  concentration  resulting  from 
chemical  reaction  and  volumetric  change  behind  the  shock  wave  is  given  by 
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where  the  positive  term  on  the  right  of  the  equality  reflects  the  sum  of  the 
chemical  rates  leading  to  species  production  and  the  negative  term,  the  rates 
for  removal.  Events  measured  in  laboratory  time  are  related  to  the  particle 
time  in  tne  shock  tube  by 
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The  above  set  of  differencial  equations  is  solved  on  an  IBM  360  computer,  using 
a  fourth-order  predictor-corrector  integration  scheme.  The  predictor-corrector 
method  was  devised  to  vary  Its  own  integration  interval  depending  on  the  esti¬ 
mated  truncation  error  available  at  each  step.  A  Runge-Kutta  method  was  used 
to  obtain  the  starting  values. (20)  However,  a  computational  difficulty  imposes 
itself  on  the  formulation  of  this  type  of  problem,  since  the  time  derivative  of 
the  area  goes  to  infinity  when  x  approaches  zero.  Teare*-®  has  overcome  this 
singularity  by  assuming  for  small  x*s  that  the  area  variation  is  a  linear  func¬ 
tion  dependent  on  the  translational/rotational  relaxation  thickness  of  the  shock 
wave.  In  the  present  program,  the  computational  difficulty  is  avoided  by  simply 
introducing  an  arbitrary  parameter  to  start  the  area  change  at  an  infinitesimal 
distance  behind  the  shock  front.  Calculations  show  that  the  results  are  vir¬ 
tually  unaffected  by  this  parameter  for  values  less  than  0.0001  of  the  chemical 
relaxation  distance. 

The  computer  program  is  general  in  that  the  Inputs  to  the  calculation  are  the 
shock  velocity,  the  state  of  the  gas  ahead  of  the  shock  wave,  and  the  chemical 
species  and  reactions  to  be  considered.  The  thermodynamic  properties  of  the 
species  used  in  the  calculation  are  automatically  furnished  by  a  computer-tape 
litrary  containing  information  on  almost  every  molecule  which  may  be  of  interest 
for  shock  tube  study.  The  computer  program  is  capable  of  solving  chemistry 
problems  with  up  to  100  reactions  and  200  species. 
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